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Abstract 

We present the order contributions to the coefficient functions corresponding to 
the longitudinal fragmentation function Fl{x,Q'^). A comparison with the leading 
order as result for Fl{x, Q^) shows that the corrections are large and vary from 44% 
to 67% in the region 0.01 < x < 0.9 at = M|. Our calculations also reveal that 
the ratio of the longitudinal and total cross section ai/atot amounts to 0.054. This 
number is very close to the most recent value obtained by the OPAL collaboration 
which obtained 0.057 ± 0.005. 



In this paper we present the 0{a'^) QCD corrections to the longitudinal fragmen- 
tation function measured in the process 

e+e- ^7,Z^i/+ "X", (1) 

where "X" denotes any inclusive final hadronic state and H represents either a specific 
charged outgoing hadron or a sum over all charged hadron species. This process has 
been studied over a wide range of energies at many different e"^ e~ colliders. For the 
most recent experimental results we refer to [1, 2]. Following the notations in [3] the 
unpolarized differential cross section of process (1) is given by 

(Pa 3 2n\ 3 . 2 c?t^L ,3 daA 

1 + cos 0) — h-sm 0— h-cosfc*— ; — . 2 



dxd cos 9 8 dx 4 dx 4 dx 

The Bj0rken scaling variable x is defined by 

Q 



x^-^ , 0<a;<l, (3) 



where p and q {q^ — Q"^ > Qi) are the four-momenta of the produced particle H and 
the virtual vector boson (7, Z) respectively. The variable 9 denotes the angle of 
emission of particle H with respect to the electron beam direction in the CM frame. 
The transverse, longitudinal and asymmetric cross sections in (2) arc defined by gt, 
aL and cx^ respectively. The latter only shows up if the intermediate vector boson 
is given by the Z-boson and is absent in purely electromagnetic annihilation. From 
the cross sections ak {k = T,L,A) one infers the fragmentation functions Fk{x,Q^) 
which are the analogues of the deep inelastic structure functions measured in deep 
inelastic lepton-hadron scattering where q is spacelike. The former are defined by 

^-a^:'mF,ix,Q^ (4) 

cr^^\Q'^) stands for the pointlike annihilation cross section presented in eq. 2.12 of [3]. 
In the framework of the QCD improved parton model the fragmentation functions 
Fk{x,Q^) can be expressed as follows 

'1 dz 



z 



D, g,^^) {z,^)+D, (|,;.^) (.,^)] ,(5)| 



where Dq and Dg are the light quark and gluon fragmentation functions as defined in 
[3] where one has summed over all charged hadron species. Furthermore the timelike 
light quark and gluon coefficient functions are given by 'Ck,q and 'Ck,g respectively 
where ji denotes the factorization scale which for convenience has been set equal to 
the renormalization scale. Notice that one can also include in (5) the contributions 
coming from the heavy quark coefficient functions and higher twist effects as has been 
done in [3]. 

The above coefficient functions have been calculated up to first order in the strong 
coupling constant as{n^) in [4, 5] where they are presented in the MS-scheme. Us- 
ing the next-to-leading order timelike DGLAP splitting functions, which have been 
computed in [6], one has made a complete next to leading order (NLO) analysis for 
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Ft{x,Q'^) in [3]. This analysis has been performed in the annihilation scheme in [3] 
where one has also included the contributions due to the charm and bottom quarks 
and higher twist effects. Here the annihilation scheme (AS) is defined by requiring 
that the total fragmentation function 

F{x, ^ Ft{x, Q^) + Fl{x, Q^), (6) 

does not recieve QCD corrections via the coefficient functions provided //^ = Q^. 
However the NLO analysis of the longitudinal fragmentation function Fl{x,Q^) is 
not complete yet because the contributions to its coefficient functions are 

scheme independent so that the scheme dependence of the two-loop DGLAP splitting 
functions is not cancelled up to 0{al). We want to fill in this gap by including the 
0{a'^) contributions to the longitudinal coefficient functions CL,i{z, Q^/yU^) {i = q, g) 
so that one can study the effects of the latter on the analysis of Fl{x, Q^) in [3]. 
The computation of ^ also leads to the 0(0^) contributions to the transverse and 
longitudinal total cross sections defined by 

from which one can derive the total cross section 

atot(g') = aT(Q') + ai(g'), (8) 

and the ratios Ufc/o'tot- Before the calculations presented below these ratios are only 
known up to 0{as) and one obtains 



Choosing as{Mz) = 0.126 [2] the quantities ctt /c^tot and Ci/crtot become 0.960 and 
0.040 respectively. These numbers have to be compared with the most recent ex- 
perimental values obtained in [1] which yield ar/cytot = 0.943 ± 0.005 and olIoxoi = 
0.057 ± 0.005. In particular the experimental value of OL/crtot is far above the theo- 
retical 0{as) prediction which indicates that higher order QCD corrections are im- 
portant. 

The calculation of the 0{a1) longitudinal coefficient functions proceeds in exactly 
the same way as done for the Drell-Yan process in [9] and deep inelastic lepton 
hadron scattering in [11]. First one computes the parton fragmentation functions 
^L,i {i — Q, g) corresponding to the process 

y ^ "a" +ai + a2 + ... + a„, (10) 

where V = 7, Z, "a" denotes the detected parton and a.i {i = 1 . . .n) stand for the 
partons of which the momenta are integrated over so that the process is inclusive 
with respect to the a^. In zeroth order of as we have the Born reaction 

V ^ V + (11) 
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which contributes to Tr.q and TA,q but not to Th.q- In next-to-leading (NLO) order 
one obtains the one-loop virtual corrections to reaction (11) and the parton subpro- 
cesses: 

V ^ + (12) 

V ^ V+g + g- (13) 

Besides J^r.i, ^A,i (^ = (1,9), the longitudinal fragmentation function Tl^x now also 
receives contributions from the above subprocesses. 

After mass factorization the coUinear divergences are removed and the objects 
turn into the coefficient function which are presented in [2], [4] and [5]. The determi- 
nation of the 0(cP'^ corrections involves the computation of the two- loop corrections 
to (11) and the one-loop corrections to (12), (13). Furthermore one has to include 



the following subprocesses: 

V ^ V+g + ^ + ^, (14) 

V ^ v'+g + g + ^, (15) 

V ''q''+q + q + q. (16) 



In reaction (16) the two anti-quarks can be identical as well as non identical. Notice 
that in the above reactions the detected quark can be replaced by the detected anti- 
quark so that in reaction (16) one can also distinguish between final states containing 
identical quarks and non identical quarks. In this paper we are only interested in the 
computation of J-L,i- Since the latter does not depend on the intermediate vector 
boson we perform our calculation for V = ^f. Notice that this also holds for J^T,i- 
After mass factorization and renormalization for which we have chosen the MS-scheme 
(indicated by a bar on the coefficient functions) the longitudinal coefficient functions 
read as follows. For the non-singlet part we have up to 0{a1) 
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where the colour factors are given by Cp = [N'^ — 1)/2N, Ca = N and Tf = 1/2 [N — 
3 in QCD). The quantity Uf denotes the number of hght flavours. The definitions 
for the polylogarithmic functions Li2(x) and Sn,p{x) can be found in [12]. In NLO 
the terms proprtional to Cp, Ca Cp and CpTf receive contributions from identical as 
well as non identical (anti) quarks in reaction (16) whereas the term proportional to 
Cp—CA Cf/2 can be only attributed to identical (anti) quarks. The singlet coefficient 
function can be written as 



^5 



(18) 



The pure singlet part denoted by Cf '^ originates from process (16) where one gluon 
is exchanged in the t-channel. It reads as follows 

cZi^, QVl^') = nf C^T, [ I 161na; + - 16a; + fx' | In ^ 

(32 16 \ 
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(19) 



Finally we have the gluonic coefficient function which is due to the parton subpro- 
cesses (13), (15). It becomes equal to 
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Besides in the MS-scheme the mass factorization can also be performed in the AS- 
scheme (see above (6)). Using the transformation formulae in [3] we obtain 



c^J(^,W) = c2,;(a:,g7/^') 
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Notice that in the above formulae we did not distinguish between the renormahzation 
and factorization scale. The distinction can be easily made if one substitutes in the 
above coefficient functions 



47r 



(24) 



where fi and fiR denote the factorization and renormahzation scale respectively. In- 
spection of the above coefficient functions reveals that there are logarithms of the 
type In x/x [k = 0, 1,2) present in the coefficient functions Cj^ ^ (19) and C^ ^ (20) 
which become large when a; — > 0. The relevance of these terms for the NLO correc- 
tions to Fl{x,Q'^) (5) will be discussed at the end of the paper. The ratios (Jk/crtot 
{k = T,L) are given by 
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(25) 



where Re+e- is defined by Re+e- — '^tot/'^^^^- Furthermore we have used in the 
derivation of (25) the equality Ck,q = Ck,q- The ratio Re+e- = o'tot/c^°^ has been 
calculated up to 0{a^g) in the literature [15]. Up to 0{a1) it is given by 



fl,, . _ = 1 + ^ CW3] + ( ^) ' [ c| { - ? } + { 



+ 
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+ nfCFTf\ 41n% + 16C(3) -22 



llln^-44C(3) 
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Substituting the coefficient functions (17)-(20) or (21)-(23) into (25) and expanding 
R'l _ (26) up to 0{al) we obtain 



0"L _ 

2023 



S)-^i^i-(S) 



33 
~2 



Cp {——)■ + CaCf 



111 24 



+ 



30 



74 



(Tt 



= 1 



(27) 
^28) 

Notice that the above perturbation series are factorization-scheme independent. They 
do however depend on the renormahzation scheme (here MS) which is indicated by /ir. 
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Since crtot = cr^*^^ Re+ e- we infer that in zeroth order, (Jtot only receives contributions 
via (7t whereas in first order of ag, (Jtot is determined by ctl. In second order, ctl as 
well as (Jt contribute to citot- 

Choosing as{Mz) = 0.126 and n/ = 5 [2] we get the following results 

— = 0.040 + 0.014 = 0.054, (29) 

— = 0.946, (30) 

which are very close to the values obtained by the OPAL-experiment [1] mentioned 
below (9). In (29) we have explicitely shown the 0{as) and 0{al) contribution 
which are represented by the numbers 0.040 and 0.014 respectively. From the latter 
we conclude that the NLO correction to oi^joiox. amounts to about 30% of the LO 
one. The agreement between the NLO theoretical predictions and the experimental 
results show that equations (27) and (28) are very suitable to determine the strong 
coupling constant provided we also include the contributions due to heavy quarks 
(here c and 6) and higher twist effects (see [3]). 

Finally we want to study the effect of the 0{p?^ corrections to the longitudinal frag- 
mentation function Fi^{x^Q'^\ Adopting the AS-scheme (see [6] above) and choosing 
the sum over the light quark fragmentation functions Dt = Dq and the gluon 
fragmentation function Dg in table 1 of [3] we have computed Fi{x^Q^) in [5] up to 
NLO. Further we took a,(Mz) = 0.126 (MS) and n/ = 5 [2]. In fig. 1 we have plotted 
the LO as well the NLO result for 0.01 < x < 0.9 and have shown for comparison 
F{x,Q'^) in (6). Notice that the data of the experiments in [1, 2] are taken in the 
range 0.01 < x < 0.9. To exhibit the NLO corrections more clearly wc have also 
plotted the ratio K = fnlo/^lo -^^ 2. The last figure shows that the O(a^) 
corrections to Fl{x,Q'^) are large in spite of the fact that as{Mz) is small. They 
amount to 45% ai x — 0.01 and increase to 67% around x — 0.1. Above x — 0.1 
the corrections decrease to 44% near x — 0.9. We have also studied the corrections 
in the region 10~^ < a; < 10~^. They are smaller than in the range 0.01 < x < 0.9 
with a minimum of 27% at a; = 7 ■ 10^^. A study of the various parts contributing 
to Fl^XjQ"^) reveals that the gluonic part in (5) dominates Fl^x^Q"^) for x < 0.2 
whereas the hght quark contribution becomes dominant for x > 0.2. This holds for 
the 0{as) as well as 0{a:^) contributions. Finally we want to make a comment on 
the contribution of the logarithmic term ln*^a;/a; mentioned under (24). They do not 
dominate the corrections when convoluted with the gluon and quark fragmentation 
functions Dg and D^. This becomes clear if one looks at the coefficient of the In^ x/x 
term appearing in the CaCf part of Ci^g in (20). Although this coefficient is neg- 
ative the actual total 0{al) contribution to F^^°(x, Q^) is positive over the whole 
x-range. 
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1 Figure captions 

Fig. 1 The leading order (LO) and next-to-leading order (NLO) corrected lon- 
gitudinal fragmentation function Fi^{x,Q'^). Dashed line: F^'^\x,Q'^)] solid line: 
fI''''°\x,Q^); dotted line : F(NLO)(x,g2) (6). 

Fig. 2 The ratio K{x, Q^) = F^^'^\x, Q^)/F^'^\x, Q^). 
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